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PREFACE 


This  report  contains  three  recent  papers  on  the  stability  £md  properties  of 
compos itionally  modulated  films. 
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INTEROIFFUSION  AND  STABILITY  Of  COHPOSITIONALLY  MODULATED  FILMS 


Frans  SPAEPEN* 

Division  of  Applied  Sciences,  Harvard  University,  Cambridge,  MA  02138 


ABSTRACT 

The  thermodynamics  and  kinetics  governing  the  stability  of  compositio- 
nally  modulated  materials  are  reviewed.  Typical  results  of  interdiffusion 
experiments  on  fee  and  amorphous  metals  are  discussed  to  demonstrate  the 
ser,si,ivity  and  special  features  of  the  technique.  Two  coarsening  irechu- 
nisms,  and  the  stability  of  orthogonal  waves  in  phase  separating  systems 
are  analyzed. 


INTRODUCTION 

Compositional ly  modulated  materials  are  produced  by  multiple  alternate 
deposition  of  very  thin  layers  of  different  composition.  The  deposition 
aethods  include  vapor  deposition  [1,2],  molecular  beam  epitaxy  [3],  sputte¬ 
ring  [4,5;,  and  chemical  vapor  deposition  [6].  If  a  structure  has  a  reaso¬ 
nably  low  free  energy  over  an  extended  continuous  composition  range,  and 
if  there  is  sufficient  intermixing  during  the  deposition  process,  the 
..ultiple  layering  process  can  result  in  the  formation  of  a  compositional  ly 
modulated  material,  i.e,  one  in  which  the  composition  profile  is  closely 
approximated  by  a  simple  harmonic  function  : 

c  •  Cjj  ♦  A  cos  ax  (1 ) 

'/here  Cg  is  the  average  composition,  x  the  distance  coordinate  normal  to 
•jie  film  plane,  A  the  modulation  amplitude  and  a  the  modulation  wavenumoer 
,3  =  2»fi;  k  is  the  modulation  wavelength).  Of  interest  here  are  materials 
with  *  between  0.2  and  10  nm.  This  composition  modulation  can  be  measured 
by  X-ray  diffraction  17, 8],  It  gives  rise  to  a  peak  at  k  =  ax,  with  intensi¬ 
ty  proportional  to  f?  .  This  peak  can  be  considered  as  a  satellite  of  tne 
(000)  forward  scattering  peak,  and  is  independent  of  the  atomic  scale  struc¬ 
ture  (crystalline  or  amorphous)  of  the  material.  In  crystal  1 ine  materials, 
satellites  also  appear  at  khid  *  a*  around  the  higher  order  Bragg  peaks 
(hkl)',  their  intensity,  however,  is  not  just  dependent  on  a’ ,  but  also  on 
the  variation  in  lattice  spacing.  Deviations  from  a  simple  harmonic  compo¬ 
sition  modulation  (eq.  (11)  can  be  detected  from  the  higher  order  peaks  at 
k  =  nax,  (n  =  2,  3,  . . . ) . 

Crystalline  compositionally  modulated  metallic  films  have  been  made 
most  successfully  in  binary  fee  systems  with  full  solid  solubility.  For 
example  :  Ag-Cu  [2],  Cu-Pd  [9],  Au-Ni  llOj.  Au-Cu  [11]  and  Cu-Ni  [12,13]. 
Usually  a  mica  substrate  was  used,  which  produced  films  with  a  strong  (ill) 
texture  that  were  fully  coherent  for  small  enough  wavelength. 

Amorphous  compositionally  modulated  metallic  films  have  so  far  only 
been  reported  for  metal -metalloid  alloys,  such  as  tPUBsSii^)  -  (FeosBij) 
]5,  14-16].  Modulated  binary  amorphous  systems  have  only  very  recently  been 
obtained;  for  example  Cu-2r  [17-18]  and  Ni-Zr  [18].  These  amorphous  modu¬ 
lated  films  differ  in  important  ways  from  the  crystalline  ones  : 
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<1)  They  can  be  produced  at  any  value  of  x ,  without  Interference  between 
the  periodicity  of  the  atonic  structure,  as  Is  the  case  at  short  modulation 
wavelengths  in  crystals.  Because  the  spatial  distribution  of  atom  centers 
In  an  amorphous  system  Is  continuous  and  the  c^osltlon  c  In  equation  (1) 
Is  the  average  one  over  the  yt-plane  at  position  x,  x  can.  In  principle, 
be  made  arbitrarily  small.  In  an  X-ray  diffractometer  experiment  with  the 
scattering  vector  k  along  x,  the  same  averaging  is  done,  and  therefore  a 
peak  at  an  arbitrarily  large  value  of  6  could  be  observed.  On  a  scale, 

however,  the  composition  gradient  can  not  be  much  greater  than  1/a  la  : 
average  interatomic  distance),  even  If  dc/dx  in  the  averaged  equation  ()) 
is  very  large  at  small  x.  This  Is  important  in  view  of  the  local  composi¬ 
tion  gradient  contributions  to  the  free  energy,  discussed  below.  In  crys¬ 
talline  materials,  the  continuum  description  must  be  corrected  when  x 
approaches  2a;  a  similar  correction  seems  required  for  amorphous  materials. 
(11)  In  many  systems,  amorphous  alloys  can  produced  over  a  much  larger 
continuous  composition  range  than  their  crystalline  counterparts.  Producing 
a  composition  modulation  In  such  a  system  may  therefore  only  be  possible 
in  the  amorphous  state. 

(iii)  Coherency  strains,  which  can  contribute  considerably  to  the  free  ene'- 

?y  of  crystalline  modulated  films  and  must  be  taken  into  account  in  the  ana- 
ysis  of  interdiffusion  [9],  are  absent,  presumably  at  least  initially,  in 
amorphous  metallic  modulated  films.  They  may  arise  from  chemical  changes 
during  the  interdiffusion  process  (19). 

(ivl  Short  circuit  diffusion  paths,  especially  grain  boundaries  (see  below), 
are  not  found  in  the  amorphous  films.  This  simplifies  the  analysis  of  the 
early  stages  of  the  interdiffusion,  and,  from  a  practical  point  of  view, 
enhances  tne  stability  and  useful  life  of  the  films. 

Artificial  compositional ly  modulated  materials  described  above  are 
thermodynamically  unstable.  They  have  an  unusually  high  density  of  inter¬ 
faces,  which  can  be  eliminated  In  two  ways  ;  homogenization  by  interdiffu¬ 
sion,  or  sharpening  of  the  composition  profile  followed  by  coarsening.  The 
thermodynamic  conditions  determining  the  evolution  of  a  modulated  structure 
will  be  discussed  In  Section  2.  Monitoring  the  Intensity  of  the  modulation 
satellites  as  a  fu'.ction  of  annealing  time  Is  the  most  sensitive  method  for 
measuring  interdiffusivity,  which  Is  especially  useful  for  the  study  of 
amorphous  metals;  this  will  be  discussed  In  Section  3.  Some  mechanisms  of 
coarsening  are  identified  and  analyzed  in  Section  4.  A  more  detailed  survey 
of  some  of  the  theory  and  experiments  discussed  here  can  be  found  in  a 
forthcoming  review  paper  (20!. 


thermodynamics 

The  theory  of  the  thermodynamics  and  stability  of  systems,  such  as  the 
modulated  films,  that  are  Inhomogeneous  on  a  very  fine  scale  was  first  fully 
formulated  by  Cahn  and  Hilliard  (21-23).  A  detailed  exposition  can  be  found 
in  the  excellent  review  paper  by  Hilliard  (24).  A  review  of  the  historical 
development  was  given  by  Cahn  In  his  1967  Institute  of  Metals  Lecture  (25  . 
In  this  section,  the  most  Important  results  are  outlined. 

In  an  Inhomogeneous  system,  the  free  energy  of  a  volume  element  is  not 
only  dependent  on  the  local  composition,  but  also  on  that  of  the  neighboring 
volume  elements  through  the  formation  of  Interfaces  with  a  composition  gra¬ 
dient.  For  an  Inhomogeneity  in  one  dimension,  x,  the  total  free  energy  of 
t.*!:  system  can  then  be  written  as  : 

F  »  ot  /(flc)  ♦  **11*^1  ***  (7) 


3 


ct  Is  the  area  perpendicular  to  ii;  f(c)  Is  the  free  energy  per  unit  volume 
of  a  homogeneous  system  of  composition  c;  k  is  the  gradient  energy  coeffi¬ 
cient.  If  F  is  minimized,  under  the  constraint  of  conservation  of  mass  the 
following  stability  conoition  is  obtained  : 


f  ■ 


3^C 

17 


13) 


where  i  is  a  constant  if  the  system  is  in  equilibrium.  Note  that  f  = 
-  Ligl/V  =  constant  (u^  :  chemical  potential  of  element  i,  V  :  molar  vo¬ 
lume)  is  the  equilibrium  condition  for  a  system  with  negligible  gradient 
energy.  If  such  a  system  is  not  in  equilibrium,  the  gradient  of  f  becomes 
the  driving  force  for  the  interdiffusion  flux,  D,  towards  equilibrium  ; 


where  M  is  the  mobility.  In  the  presence  of  gradient  energy  contributions, 
the  driving  potential  for  interdiffusion  becomes  a  in  equation  (3)  : 


J 


3X 


15) 


Combination  of  equations  (3)  and  (5)  with  the  conservation  of  mass,  3c/3t= 
-  V  .  div  0,  gives  the  diffusion  equation  : 


•t 


3^C 
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16) 


where  M.  f  and  <  are  assumed  constant  to  obtain  a  linear  equation.  The 
bulk  interdiffusion  coefficient  is  : 

13  =  M  f"  V  17) 


The  evolution  of  a  composition  modulation  wave  is  obtained  by  inserting 
equation  11)  in  the  diffusion  equation,  which  yields  : 


A 


exp  1  -  Ob 


2 


18) 


'he  relative  rate  of  change  in  amplitude  is  defined  as  an  amplification 
factor  : 


3A 

5t 


19) 


where 


is  the  wavelength-dependent  interdiffusion  coefficient. 


The  growth  or  decay  of  a  modulation  wave  now  depends  on  the  sign  of 
R.  Three  different  cases,  illustrated  on  Figure  1,  can  be  distinguished. 
To  facilitate  the  interpretation  it  is  useful  to  consider  the  predictions 
of  the  regular  solution  model  for  f"  and  ».  For  an  equiatomic  binary 
mixture,  with  a  molar  energy  of  mixture  au  ; 


f"  ,  1  (  RT  -  2aU  . 

V 

The  gradient  energy  coefficient  is  (21) 


V 


111) 

(12) 
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Figure  la  (f  <  0,  «  >  0)  corresponds  to  a  phase  separating  system  (iU  >0) 
inside  the  spinqdal.  Modulations  with  wavelengths  greater  than 
=  2  «  (-  2  c/f")'  grow;  this  is  the  classical  case  of  spinodal  decomposi¬ 
tion,  corresponding  to  uphill  diffusion  (D  <  0,  eq.  (7)).  Modulations  with 
i  <  decay,  due  to  their  high  density  of  interfaces  which  make  a  positive 
gradient  energy  contribution. 

Figure  lb  (f"  >  0,  «  >  0)  corresponds  to  a  phase  separating  system  laU  >0) 

outside  the  spinodal.  The  diffusion  coefficient  is  positive  for  all  wave¬ 
lengths,  and  increases  with  decreasing  wavelength  due  to  the  increasing 
driving  force  for  homogenization  from  the  interfacial  energy.  Figure  Ic 
(f"  <  0,  r  <  0)  corresponds  to  an  ordering  system  (all  <  0),  and  is  the  exact 
negative  of  Figure  la.  At  long  wavelengths  t»  >  i^),  the  system  homogenizes 
(D  >  0,  eq.  (7)).  Modulations  with  short' wavelength  (»  <  are  predicted 
to  grow;  since  in  crystalline  systems  is  usually  on  the  order  of  the  ato¬ 
mic  spacing,  the  development  of  very  short  wavelength  “phase  separation” 
can  be  interpreted  as  an  ordering  process.  It  should  also  be  kept  in  mind 
that  the  continuum  approximation  used  above  must  break  down  at  small  k. 
Corrections  for  the  discreteness  of  the  crystalline  lattice  have  been  made 
by  Cook  et  al.  (261.  The  fourth  possibility  (f"  <  0,  «  <  0)  is  an  unlikely 
one.  In  fact,  it  can  not  occur  in  the  regular  solution  approximation  (equa¬ 
tions  (11)  and  (12)). 


Since  the  lattice  parameter  of  an  alloy  depends  on  its  composition 
(n  >  d(ln  a)/dc),  matching  of  the  lattice  planes  along  a  modulation  wave 
introduces  coherency  strains.  The  free  energy  of  eq.  (2)  now  acquires  an 
additional  term  [22,23  |  ; 

F  =  o<  /  :f(c)  ♦  X  (|-|)^  (c  -  Cjjl^ldx 


(13) 


where  Y  Is  a  modulus  that  depends  on  the  crystal  synnetry  and  the  direction 
of  the  modulation  wave,  for  an  isotropic  continuum  [221  : 


Y  =  E/(l  -  '<)  (14) 

E  :  Young's  modulus;  “  :  Poisson’s  ratio. 

to  a  result,  the  interdiffusion  coefficient  now  contains  an  additional  posi¬ 
tive  (i.e.  homogenizing)  term  : 


=  6  ;i  * 


1  *  2M"^Y 


(15) 


fete  that  f"  now  must  be  more  negative  to  cause  spinodal  decomposition  than 
in  the  absence  of  coherency  strains.  The  condition  f"  <  -  2n^'<  defines  the 
"coherent  spinodal". 


INTERDIFFUSION 

2 

Since  the  intensity.  I,  of  the  (000)  satellite  is  proportional  to  A  , 
(as  are  the  (hkl)  satellites  in  the  absence  of  coherency  strains)  is  the 
slope  of  a  plot  of  In  I  vs  t,  according  to  eq.  (9),  equal  to  -  2  D»s^.  Sucn 
a  plot  usually  shows  an  initial  transient  with  a  decreasing  slope  before 
settling  down  to  a  “terminal  slope"  from  which  the  interdiffusivity  is  taken. 
Cook  and  Hilliard  [21  have  attributed  these  transients  to  increased  homoge¬ 
nization  in  moving  grain  boundaries  as  a  result  of  texture-driven  gram 
growth  I  see  Figure  2) . 

F'x.pu-te  2  : 
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In  interdiffusion  experiments  on  fcc-metallic  compositional ly  modulated 
laterials.  the  wavelength  dependence  of  IH  was  found  to  obey  equations  (10) 
or  (151,  and  the  first  direct  determinations  of  the  gradient  energy  coeffi¬ 
cient,  X  ,  were  made.  Cook  and  Hilliard  [2)  studied  the  Ag-Au  system,  in 
which  the  coherency  strain  effects  are  negligible  (n  =  0.15  t).  A  fit  of 
their  data  to  eq.  (10)  gave  x  =  -  2.6  x  10*''  J.m'’  ,  a  negative  value  as 
expected  for  an  ordering  system  (Figure  Ic).  In  the  temperature  range  200- 
2S0°C  they  measured  bulk  interdiffusivities,  D,  as  low  as  lO"^^  m‘.s”', 
which  were  in  perfect  agreeiKnt  wi th  an  extrapolation  of  high  temperature 
(i  900°C)  data  (D  >  10"'*  nr  .s'  ),  which  clearly  represent  lattice  diffu¬ 
sion.  The  modulated  film  technique  is  therefore  not  only  the  most  sensitive 
one,  but  it  also  allows  determination  of  the  lattice  diffusivity  in  a  tempe¬ 
rature  regime  (T  <  2T)i|/3t  where  other  techniques  would  only  detect  short- 
circuit  diffusion  along  grain  boundaries  or  dislocations.  The  reason  for 
this  is,  of  course,  that  the  diffusion  distance  for  lattice  diffusion  (i/2) 
IS  much  lesser  than  either  the  grain  size  or  the  distance  between  disloca¬ 
tions,  so  that  the  fraction  of  the  material  homogenized  by  short  circuit 
diffusion  is  negligible  to  that  homogenized  by  lattice  diffusion.  The  only 
snort  circuit  mechanism  that  can  homogenize  an  appreciable  amount  of  mate¬ 
rial  is  that  involving  moving  grain  boundaries,  as  discussed  above  (see  Fi¬ 
gure  2).  This  mechanism  also  explains  the  very  high  effective  diffusivities 
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observed  In  the  earliest  coapositlonally  Modulated  Materials  [lli  Mhich  had 
a  fine-grained  polycrystalline  structure. 

The  effect  of  the  coherency  strain  on  the  dlffuslvity  was  confirmed 
by  Philofsky  and  Hilliard  (9)  for  the  Cu-Pd  syste_.  (n  =  7  1).  By  comparing 
the  intensities  of  the  satellites  at  kjj)  ♦  a  and  k  .  p  -  B,  they  establisheo 
that  for  full  coherence  the  mooulation  wavelength 'had  to  be  less  than  2.E 
kw;  beyond  that,  a  gradual  loss  of  coherence  was  observed  with  full  incohe¬ 
rence  beyond  3.8  mm.  The  wavelength  dependence  of  D>  was  fpund  to  obey 
equation  (15)  for  k  <  2.8  mm;  beyond  that  a  gradual  decrease  of  Dj  with  ; 
was  observed,  corresponding  with  the  decrease  in  the  positive  strain  term 
in  equation  (15) . 


The  modulated  film  technique  has  special  advantages  in  the  study  of 
diffusion  in  metallic  glasses. 

(i)  Since  Metallic  glasses  crystallite  rapidly  if  their  dlffuslvity  is 
greater  than  about  10"”“  m^.s"',  sensitive  detection  techniques  are  required. 
The  solution  of  all  other  methods  for  Measuring  the  dlffuslvity  is  determi¬ 
ned  by  the  spatial  resolution  of  the  "slicing*  technique  used  to  analyze 
a  single  diffusion  junction.  Specialized  techniques,  such  as  sputter  pro¬ 
filing,  Rutherford  backscattering  or  nuclear  reactions  must  be  used  on  the 
metallic  glasses,  and  thev  are  IlMited  to  only  a  few  orders  of  magnitude 
of  D.  Below  10"23  m2.s"',  only  the  modulated  film  technique  is  available. 
Its  sensitivity  is  a  result  of  its  being  composed  of  several  hundred  diffu¬ 
sion  junctions,  so  that  very  small  changes  in  the  composition  profile, 
corresponding  to  a  diffusion  distance  of  no  more  than  10  pm,  can  be  detected 
by  the  X-ray  technique. 

(ii)  Since  glasses  are  thermodynamically  unstable,  they  undergo  a  continuous 
series  of  transformations  to  states  of  lower  free  energy.  This  process  of 
structural  relaxation  affects  all  physical  properties,  but  most  strongly 
the  atomic  transport  properties  such  as  viscosity  (27]  and  dlffuslvity 
[is;.  Since  the  modulated  film  technique  Is  non-destructive  it  can  be  used 
conveniently  to  measure  the  time  dependence  of  the  dlffuslvity .  Metallic 
glasses  show  a  very  large  transient  In  their  In  1  vs  t  curve,  which  can  only 
be  attributed  to  structural  relaxation  (15);  Indeed,  for  a  (Pdgj Si Ign / 
(FegsBiglgo  film,  the  quantity  (t)  was  observed  to  Increase  linearly 
with  time,  exactly  as  one  would  expect  from  a  scaling  relation  with  the 
viscosity,  which  has  also  been  observed  to  Increase  linearly  with  annealing 
time  127  I .  The  technique  also^allows  Iso-configurational  measurements  of 
the  temperature  dependence  of  6.  After  sufficient  annealing,  structuraj 
relaxation  becomes  slow  enough  so  that  upon  cycling  of  the  temperature,  D 
can  be  reproduced.  The  Iso-configurational  activation  energy  for  diffusion 
measured  this  way  was  very  close  to  that  of  the  viscosity  in  a  similar 
system  ( 15,27] . 

(iil)  Since  interdiffusion  Is  driven  by  a  chemical  driving  force,  the  modu¬ 
lated  film  technique  allows  Insight  Into  the  thermodynamics  of  metallic 
glasses,  which  are  difficult  to  obtain  otherwise.  Canmerata  and  Greer  116) 


have  studied  the  wavelength  dependence  of  Dx  In  the  tP<i85^M5 150^*^*85^15*50 
system.  Their  results  are  redrawn  on  Figure  3,  together  with  the  fit  to 

equation  (10).  The  data  points  were  taken  from  the  terminal  slopes  of 


In  1  vs  t  curves.  Since  the  films  had  identical  fabrication  and  thermal 


histories,  and  exhibited  similar  relaxation  kinetics,  the  data  were  assumed 
to  correspond  to  the  same  degree  of  relaxation.  Nevertheless,  part  of  the 
scatter  in  the  data  may  still  be  due  to  the  difficulty  in  reproducing  an 
identical  structural  state.  The  system  Is  clearly  a  homogenizing  one. 
Since  the  bulk  Interdiffusivlty.  _0,  is  positive,  f”  >  0  according  to  equa¬ 
tion  (7).  Since  the  slope  of  vs  B  Is  negative,  c  must  be  negative, 
which,  according  to  the  regular  solution  approximation  (equation  (12)), 
corresponds  to  an  ordering  system.  The  equiatomic  heat  of  mixing,  a u,  for 
Pd  and  Fe,  which  are  the  heaviest  elements  whose  Interdiffusion  is  mainly 
being  observed  by  the  X-ray  technique,  has  Indeed  been  calculated  to  be 
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-  6  kJ  .  mole"'  126],  and  ordered  phases  occur  in  the  crystalline  state. 
A  self-consistent  calculation  of  aU  from  Figure  3  and  equations  (11)  and 
(121,  however,  is  not  possible.  It  is  probably  necessary  to  extend  the  cal¬ 
culation  of  «  beyond  the  nearest  neighbor  approximation  [291. 

The  scaling  relation  between  diffusivity,  D,  and  viscosity,  n,  [301is; 


where  L  is  a  characteristic  length.  In  the  liquid  state,  and  around  the 
glass  transition  temperature,  L  =  3»a,  corresponding  to  the  Stokes-Einstein 
relation.  To  use  the  value  of  Figure  4,  it  is  necessary  to  correct  for  the 
chemical  driving  force  (241  ; 

6  =  D  (17) 

or,  with  c  =  0.5  and  the  regular  solution  approximation  of  equation  (11) 
0-0(1-  ZaU/RT).  Using  the  calculated  value  all  ■  -  6  ItJ.mole  (281,  this 
gives  :  6  =  3.8  0.  Combined  with  earlier  results  on  the  viscosity  [15,)6, 
27;,  this  gives  in  equation  (16)  ;  I  =  0.06  a.  This  is  still  150  times 
smaller  than  the  Stokes-Einstein  value,  but  closer  to  it  than  an  earlier 
estimate  without  the  driving  force  correction  [16].  It  should,  of  course, 
be  kept  in  mind  that  this  analysis  is  only  approximate,  since  there  is 
insufficient  information  to  take  possible  chemical  effects  of  the  metalloids, 
B  and  Si,  on  the  interdiffusion  of  Fe  and  Pd  into  account. 


COARSENING 

!n  a  phase  separating  system  inside  the  spinodal  (see  Figure  la),  a 
composition  modulation  of  sufficiently  large  wavelength  (k  ’  «  c)  grows. 
This  has  been  observed  in  modulated  films  of  the  Au-Ni  [101  and  Cu-Ni  [12, 
13  systems.  The  modulation  amplitude  continuous  to  increase  until  it 


*  . 


reaches  the  equillbrlua  coeiposltlon  at  the  annealinq  teaperature.  At  the 
same  time,  ‘squaring*  of  the  composition  profile  occurs,  and  continues  until 
each  point  Is  at  one  of  the  tMO  equilibrium  coe^MSItlons,  except  for  the 
Interfaces,  which  must  have  a  finite  width  due  to  the  gradient  energy  iZI]. 
This  final  state  In  the  evolution  of  a  growing  modulation  is  called  a  "sta¬ 
tionary  state".  It  Is  clear  that  the  evolution  towards  this  state  can  no 
longer  be  described  by  the  linear  diffusion  equations  (6),  since  k,  f“  and 
M  can  only  be  assumed  constant  If  the  amplitude.  A,  is  small.  For  larger 
amplitudes,  the  composition  dependence  of  f  can  be  taken  into  account  witn 
a  composition  dependent  diffusion  coefficient,  which  is  usually  done  by  a 
Taylor  expansion  in  ac  =  c  -  c^  : 

B  *  D°  +  6'  (ac)  ♦  J  B“  (ac)^  118) 

The  coefficients  are  proportional,  respectively,  to  the  second,  third  and 
fourth  derivatives  of  f,  and  are  therefore  related  to  the  spinodal  and  equi¬ 
librium  compositions  [13]  .  Analyses  based  on  the  resulting  non-1inear 
diffusion  equation  have  been  made  by  de  Fontaine  [31  I,  Cahn  [  3Z  1,  tanger 
[33]  and  Tsakalakos  [12,13,331.  In  the  non-linear  regime,  the  Fourier 
components  of  a  composition  profile  are  coupled;  as  a  result,  higher  harmo¬ 
nics,  leading  to  "squaring”  of  an  Initially  purely  sinusoidal  wave  and  to 
saturation  of  growth  at  a  stationary  state,  are  introduced.  By  fitting  the 
experimentally  observed  evolution,  as  a  function  of  annealing  time,  of  the 
fundamental  and  the  harmonics  to  the  theory,  the  coefficients  In  equation 
(18)  can  be  determined  [12,13,351. 

The  stationary  state  Is  still  thermodynamically  unstable  due  to  its 
density  of  Interfaces.  The  free  energy  can  therefore  be  lowered  further 
by  a  continuous  process  of  coarsening.  Since  this  process  is  driven  by  the 
surface  tension,  a,  it  is  useful  to  recall  Cahn  and  Hilliard's  [21]  result 
for  a  diffuse  interface  : 


a 


(19) 


where 

t  •  I  (20) 

The  width  of  the  interface  is  2t.  In  the  classical  Llfshitz-Slyorov-Hagner 
theory  of  Ostwald  ripening  the  coarsening  rate  of  a  population  of  spheri¬ 
cal  precipitates  in  the  diffusion-limited  regime  is  given  by  [361  ; 

-  C,V  m  o  t  (21) 

>d)ere  R  is  the  average  particle  size,  Rg  the  Initial  average  particle  size, 
and,  C;  a  proportionality  constant.  Langer  [33]  has  developed  a  theory  for 
Identifying  some  of  the  stationary  states  In  one-  and  three-dimensional 
spinodally  decomposing  systems,  and  has  demonstrated  this  instability.  From 
the  rate  of  decay  of  these  states,  he  was  able  to  predict  the  coarsening 
rate.  For  a  three-dimensional  system  of  spherical  particles,  he  predicts; 

R^  -  Rg  '  C2f"  V  M  t  t  (22) 

where  C2  Is  again  a  proportionality  constant.  Taking  into  account  eq.  (19), 
the  functional  dependence  on  time  and  surface  tension  In  equations  (21)  and 
(22)  are  seen  to  be  the  same. 


For  a  one-dimensional  system,  such  as  the  modulated  structures  consi- 
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dered  here,  he  hes  Identified  ts  ‘uniform*  cotrsening  process.  Illustrated 
In  figure  4,  In  which  alternating  layers  of  one  component  thicken  and  thin 
down,  leading  to  an  eventual  doubling  of  the  wavelength.  The  coarsening 
rate,  expressed  as  a  continuous  process.  Is  then  : 


■  7 


In  (I  - 


BfVHt 
“2 — 


e 


- 


(23) 


where  Jgls  the  Initial  wavelength. 


A  second  mechanism  for  coarsening  In  modulated  films  involves  the  mo¬ 
tion  of  fault  lines  (i.e.  edges  of  missing  planes),  as  observed  In  the  coar¬ 
sening  of  lamellar  eutectics.  As  Illustrated  on  Figure  5,  the  half  plane 
recedes  by  diffusion  of  atoms  from  Its  tip  to  the  thickening  adjacent  la¬ 
mellae.  A  simple  analysis  of  the  diffusional  flow  driven  by  the  difference 
in  curvature  between  the  tip  (positive)  and  the  adjacent  lamellae  (negative), 
gives  for  the  coarsening  rate  (37)  : 


1  -  X  a  C,  0-  !  M  a  t 
0  3  f 


(24) 


where  Cg  is  a  geometrical  constant,  estimated  at  about  300;  Of  Is  the  fault 
line  density  (i.e.  the  number  of  configurations  of  Figure  5  per  unit  area). 
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Although  this  fault  mechanism  Is  likely  to  dominate  the  uniform  one, 
due  to  the  higher  driving  force.  It  Is  still  Interesting  to  compare  equa¬ 
tions  (23)  and  (24).  The  linear  dependence  on  time  of  the  fault  mechanism 
is  found  only  approximately  In  the  uniform  mechanism  of  equation  (23),  and 
only  if  X  w  t,  which  Is  the  macroscopic  limit.  The  linear  proportiona¬ 
lity  of  the  coarsening  rate,  dk/dt,  with  the  surface  tension  found  In  the 
fault  mechanism,  however,  does  not  have  a  simple  equivalent  In  the  uniform 
mechanism.  From  equation  (23)  It  Is  seen  that  dx/dt  Indeed  Increases  with 
t  (proportional  to  o;  equation  (191)  In  the  macroscopic  limit  (x  »c), 
b.t  not  linearly.  For  short  wavelengths  (  x  <  t/2),  the  coarsening  rate, 
even  decreases  with  Increasing  t.  presumably  oue  to  Increasing  overlap  be¬ 
tween  the  diffuse  interfaces  of  adjacent  layers.  It  seems  Hkely  that  the 
fault  mechanism  also  needs  to  be  corrected  In  the  short  wavelength  limit 
(X  <  0. 
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In  Mdulated  materials  of  phase  separating  systems  inside  tne  spinodal, 
one  also  must  consider  the  possibility  of  additional  phase  separation  in 
directions  other  than  these  of  the  original  modulation,  as  illustrated  on 
figure  6.  This  situation  is  different  from  that  usually  encountered  in 
spinodal  decomposition,  wnere  the  two  orthogonal  waves  usually  have  similar 
amplitudes.  In  a  modulated  film,  the  pre-eaistence  of  a  large  amplitude 
modulation  affects  the  stability  of  small  waves  orthogonal  to  it.  This  will 
be  demonstrated  below  on  a  simple  example. 

Consider  a  two-dimensional  superposition,  as  in  Figure  6  : 

ac  •  c  -  Cjj  «  cos  cos  a^y  (25) 

The  total  free  energy  of  the  system  is  given  by  the  two-dimensional  equiva¬ 
lent  of  equation  (2)  : 

F  • (flc)  ♦  fUc)^)  dxdy  (26) 

The  bulk  free  energy  must  have  at  least  one  non-parabo1ic  term,  for  example: 

f(c)  =  -  2  A  (ac)^  ♦  ^  B  (ac)*  (T?) 

The  total  free  energy  per  unit  volume  is  then  : 

‘  -  5  'Cx  <1^  *  I?  *  2  (28) 

From  thiv  equation  it  is  imediately  clear  that  if  the  amplitudes  are  small 
(C*  «  c‘l,  or  if  the  free  energy  is  parabolic  (8  »  0,  linear  diffusion 
equation),  the  two  waves  are  independent  and  have  the  same  stability  crite¬ 
rion  (aF/aCt  *  0),  s2  <  A/2x.  which  is  identical  to  the  expression  derived 
from  equation  (10).  The  case  of  interest  here  is  :  B  A  0,  and  Cy  «  Cx, 
so  that  higher  powers  of  Cy  can  be  neglected.  The  stability  criteria  for  the 
two  waves  are  now  different  : 

Sx  <  (A  -  I  BC^)/2x  (29) 


Comparison  of  equations  (29)  and  (30)  shows  that  the  spectrum  of  new  small 
amplitude  orthogonal  waves.  Cy .  is  not  only  more  restricted  than  that  of 
waves  forming  from  a  homogeneous  solution,  but  even  than  that  of  the  large 
amplitude  pre-existing  modulations. 
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A  slnlltr  quAlitulve  conclusion  can  be  reached  froa  Cahn's  non-linear 
analysis  of  the  later  stages  of  spinodal  decoapositlon  [321,  In  which  he 
derives  expressions  for  the  growth  rates  of  coupled  orthogonal  waves  of  a 
particular  wavenumber,  a,  by  successive  apprcKlmatlon  In  the  problem  of 
Interest  here,  the  spectrum  of  the  first  approximation  Is  dominated  by  a 
single  wave  of  amplitude  Cx-  The  amplification  factors,  at  t  >  0,  for  waves 
in  the  x  and  y  directions  in  the  second  approximation  are  then  : 


ac' 


3  0“ 


c;  1 


o  .  >  '’‘■y  -  0  ( w  3  ^  ,2  , 
"y  -  c:  ar-  ■ " '  ’  *  2  r  * 


(31) 

(32) 


where  ^  and  C'  are  the  amplitudes  of  the  waves  in  Cahn's  second  approxima¬ 
tion,  IT  is  the  amplification  factor  In  the  linear  regime  (equation  (9)), 
bj  is  the  Interdiffusion  coefficient  of  equation  (10)  (negative),  0“  is  the 
quadratic  coefficient  in  equation  (18),  which  is  proportional  to  f'X  *  6B. 
Again,  the  growth  rate  of  the  orthogonal  wave  is  only  affected  in  the  non¬ 
linear  regime  (B  f  0),  and  if  the  original  modulation  is  large,  tn  that 
case,  its  growth  rate  is  not  only  less  than  in  the  linear  regime,  but  also 
less  than  that  of  the  original  modulation. 


CONCLUSIONS 

Artificial  compositional ly  modulated  materials  are  inherently  unstable. 
The  modulation  disappears,  either  by  homogenization  (in  long  wavelength  or¬ 
dering  systems,  in  phase  separating  systems  outside  the  spinodal  or  for 
short  wavelengths  inside  the  spinodal),  or  by  initial  sharpening  of  the  pro¬ 
file  to  a  stationary  state  (in  long  wavelength  phase  separating  systems  or 
short  wavelength  ordering  systems)  followed  by  coarsening. 

Measurements  of  the  rate  of  homogenization  are  a  very  sensitive  method 
for  determining  interdiffusion  coefficients.  For  crystalline  materials, 
the  method  allows  determination  of  the  lattice  diffusion  coefficient  at  very 
low  homologous  temperatures.  For  metallic  glasses,  it  permits  time  depen¬ 
dent  measurements  of  D,  which  are  essential  to  take  into  account  structural 
relaxation. 

Two  coarsening  mechanisms,  a  uniform  one,  and  one  governed  by  the  mo¬ 
tion  of  fault  lines,  have  been  discussed.  Further  investigation  into  the 
coarsening  rate  of  the  fault  mechanism  for  wavelengths  on  the  order  of  the 
interface  thickness  seems  desirable.  The  presence  of  an  initial,  large 
amplitude  modulation  wave  stabilizes  the  structure,  through  the  non-linear 
effects  in  the  diffusion  equation,  against  the  growth  of  certain  waves 
orthogonal  to  the  original  one.  A  full  analytic  or  numerical  treatment  of 
this  question,  especially  taking  into  account  crystal  anisotropy,  would  be 
of  interest. 
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Interdiffusion  in  Si/Ge  amorphous  multilayer  films 

S.  M.  Prokes  and  F.  Spaepen 

Division  o/ Applied  Sciences.  Harvard  University.  29  Oxford  Street.  Cambridge.  Massachusetts  02138 
(Received  12  April  1985;  accepted  for  publication  10  May  1985) 

Multilayered  or  compositionally  modulated  amorphous  Si/amorphous  Ge  films  with  a  repeat 
length  of  5.83  nm  have  been  fabricated  using  ion  beam  sputtering.  The  interdiffusion  coefficient 
Di  was  determined  by  measuring  the  intensity  of  the  x-ray  satellite  arising  from  the  modulation  as 
a  function  of  annealing  time.  The  interdiffusion  was  found  to  be  relatively  rapid,  in  that  it  could  be 
measured  easily  without  crystallization  occurring.  The  temperature  dependence  in  the  range  550- 
630  K  is  described  by  =  1.07  X  10' exp  (  -  1.6  eVAT)  m's” '. 


Atomic  diffusion  in  amorphous  semiconductors  has  not 
been  studied  until  very  recently,  and  the  only  data  available 
so  far  are  for  impurity  diffusion. '  The  diffusion  of  the  cova¬ 
lent  random  network  formers  themselves,  i.e.,  the  Si  or  Ge 
atoms,  has  so  far  not  been  investigated.  The  main  problem 
here  is  the  thermal  stability  of  the  amorphous  phase.  Above 
about  700  K  for  amorphous  Ge^'^  and  900  K  for  amorphous 
Si,^  rapid  crystallization  sets  in.  It  is  therefore  a  priori  not  at 
all  clear  that  the  diffusion  rate  of  the  network  atoms  should 
ever  be  fast  enough  to  be  measurable. 

The  most  sensitive  technique  available  for  measuring 
diffusivities  makes  use  of  multilayered  or  compositionally 
modulated  films  with  a  repeat  length  of  a  few  nanometers. 
By  monitoring  the  intensity  of  the  x-ray  reflection  resulting 
from  the  modulation  as  a  function  of  annealing  time,  inter- 
diffusivities  as  tow  as  10~^^  m^$~'  can  be  measured.  This 
technique,  originally  developed  for  crystalline  materials,'*^ 
has  been  used  successfully  to  measure  diffusivities  in  amor¬ 
phous  metals.^*’  Since  it  is  nondestructive,  the  technique 
allows  successive  measurements  on  the  same  sample.  This  is 
particularly  important  in  amorphous  materials,  which  ex¬ 
hibit  structural  relaxation  that  affects  all  physical  proper¬ 
ties,  and,  in  particular,  the  atomic  transport  coefficients.^*  '*’ 
The  compositionally  modulated  film  technique  can  there¬ 
fore  be  u^  to  monitor  the  extent  of  structural  relaxation, 
and  to  ascertain  that  measurements  as  a  function  of  tempera¬ 
ture  are  carried  out  under  isoconfigurational  conditions  (i.e., 
with  negligible  relaxation). 

Multilayered  films  of  (hydrogenated)  amorphous  Si  (a- 
Si)  and  amorphous  Ge  (n-Ge)  have  been  made  before,"  but 
so  far  they  have  not  been  used  for  the  study  of  interdiffusion. 

The  u-Si/o-Ge  compositionally  modulated  films  used  in 
the  present  study  were  prepared  using  an  ion  beam  sputter¬ 
ing  system  consisting  of  an  ion  gun  Uut  sputters  from  alter¬ 
nating,  motor  driven  targets.  The  thickness  of  the  layers  is 
controUed  by  a  monitoring  crystal.  The  1-^m-thick  films 
were  sputtered  onto  a  copper  substrate  and  had  an  average 
composition  of  Si^Ge^,  as  determined  by  electron  micro- 
probe  analysis.  The  layer  repeat  length  (or  compositional 
wavelength),  determined  by  x-ray  diffraction,  was  5.83  nm. 
Figure  1  shows  x-ray  traces  of  the  first,  second,  and  third 
order  satellites  about  the  (000)  reflection.  The  measured  full 
width  at  half-maximum  for  the  first-order  satellite  was 
d  (261)  =  0.09*  for  Cri<r„  radiation.  This  corresponds  to  a  7% 
variation  in  the  layer  thickness  during  deposition.  The  sam¬ 
ples  were  examined  by  transmission  electron  microscopy. 


and,  from  the  diffraction  pattern  and  dark-field  images,  were 
found  to  be  completely  amorphous.  Differential  scanning 
calorimetry  (DSC)  measurements  were  also  performed  to 
check  the  thermal  stability.  At  a  scan  rate  of  40  K/min,  two 
exothermic  peaks  were  observed.  The  first,  at  800  K,  corre¬ 
sponds  to  the  amorphous-to-crystalline  transition  of  Ge;  the 
second,  starting  at  880  K,  can  be  attributed  to  the  crystalliza¬ 
tion  of  a-Si.  These  crystallization  temperatures  of  a-Gt  and 
a-Si  are  consistent  with  the  results  of  Chen  and  Turnbull^ 
and  Donovan  tt  al.^ 

One  film  was  annealed  in  a  furnace  that  was  first  evacu¬ 
ated  and  then  backfilled  with  helium.  The  temperature  was 
monitored  using  a  thermocouple  connected  to  a  chart  re¬ 
corder,  so  that  corrections  for  the  heat-up  time  could  be 
made.  After  each  anneal,  the  decay  of  the  x-ray  satellites  was 
determined  on  a  C.E  horizontal  diffractometer,  using  CtK„ 
ladiation.  X-ray  scans  were  also  made  at  appropriate  values 
of  26  to  detect  possible  crystallization  during  the  anneals. 
No  crystalline  Bragg  peaks  were  observed  during  the  entire 
experiment.  Figure  2  shows  the  decay  of  the  integrated  in¬ 
tensity  of  the  first-order  satellite  as  a  function  of  armealing 
time  for  different  temperatures.  The  decay  of  the  intensity,  I, 
is  related  to  the  interdiffusion  coefficient  by^ 


where  A  is  the  composition  modulation  wavelength  and  is 
the  initial  intensity. 

The  initial  aimeal  was  performed  at  630  K,  and  further 
measurements  were  made  at  progressively  lower  tempera- 
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FIG.  2.  Decay  of  Ihe  intensily  of  the  first  satellite  as  a  function  of  annealing 
lime  and  temperature. 

tures.  As  discussed  by  Greer  el  al.,*  this  procedure,  com¬ 
bined  with  a  sufficiently  long  first  anneal,  establishes  isocon- 
figurational  conditions.  A  plot  of  the  logarithm  of 
normalized  intensity  versus  time  for  each  anneal  is  shown  in 
Fig.  2.  The  interdiffusivity  at  630  K  was  measured  again 
following  the  series  of  experiments  at  lower  temperatures, 
and  was  found  to  be  only  3%  tower  than  the  first  value.  Such 
reproducibility  demonstrates  that  the  measurements  are,  in 
fact,  isoconfigurational. 

A  small  amount  of  nonlinearity,  most  likely  caused  by 
structural  relaxation,"  was  seen  during  the  early  part  of  the 
first  anneal  at  630  K,  after  which  the  decay  became  linear. 
Tlie  interdiiTusivity  was  determined  for  each  tempera¬ 
ture  from  the  slopes  of  the  lines  on  Fig.  2  according  to  Eq.  1 1 ). 

The  results  are  listed  in  Table  1.  As  illustrated  by  Fig.  3, 
the  temperature  dependence  can  be  described  by  Arrhenius 
law  with  an  activation  energy  of  l.S  X  lO’  J  mole~ '  and  a 
pre-exponential  factor  Do=  1.07  X  10“'®  m*s“'.  An  ex¬ 
trapolation  of  the  data  to  700  K  yields  D;^  =  6.8  X 10“^^ 
m^s“',  in  good  agreement  with  the  range  10“^'-10“” 
m^  s~  ’,  estimated  by  Tsuei  from  the  observation,  by  means 
of  the  resistivity  change,  of  the  crysullization  of  an  a-Si/n- 
Ge  multilayer. 

The  interdiffusivities  measured  in  the  present  work  are 
remarkably  high,  which  is,  of  course,  the  reason  why  they 
could  be  measured  without  any  crystallization  occurring. 
This  is  also  apparent  from  a  comparison  with  diffusion  in 
crystals.  At  630  K,  the  interdiffusivity  in  the  amorphous 
phase  is  3.86  X  10“”  m^s"',  whereas  the  only  relevant 
crystalline  dau  available  give  a  diffusivity  of  Ge  in  crystal¬ 
line  Si  or  Ge,  extrapolated  to  this  temperature,  of  only 
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FIG.  3.  Temperature  dependence  of  the  interdilTusiviiy.  D, 

4.35  X  10“”  m‘s"'  and  5.3  X  10”"'  m's”',  respective¬ 
ly.'"''’  The  high  interdiffusivity  measured  in  the  amorphous 
phase  is  certainly  not  a  result  of  short-circuit  diffusion  along 
grain  boundaries  or  dislocations.  Even  in  crystals,  it  is  the 
bulk  diffusivity  that  governs  the  interdiffusion  in  modulated 
films,  due  to  the  modulation  repeat  length  being  much 
shorter  than  the  distance  between  lattice  defects.  The  only 
effective  short-circuit  paths  identified  in  crystalline  modu¬ 
lated  films  are  moving  grain  boundaries,’  and  they  are  clear¬ 
ly  not  found  in  amorphous  materials. 

The  microscopic  mechanism  for  the  diffusion  of  the 
network  atoms  is  as  yet  unclear.  The  low  pre-exponential 
could  result  from  the  diffusion  being  governed  by  pre-exist¬ 
ing,  nonequilibrium  defects,  such  as  broken  bonds,  that  are 
relaxing  out  very  slowly.  Since  a  normal  pre-exponential  for 
substitutional  diffusion  would  be  aboi  I0""-10"*  m*  s " 
a  defect  concentration  of  about  10“’  pe,  atom  would  be  re¬ 
quired.  The  activation  energy  is  close  to  the  Ge-Ge  bond 
energy;  it  is  conceivable  that  a  diffusive  jump  of  Si  atoms  in 
Ge,  which  presumably  dominates  the  interdiffusion  process, 
requires  the  breaking  of  an  extra  matrix  bond  next  to  the  pre¬ 
existing  defect. 

Before  the  diffusion  mechanism  can  be  positively  identi¬ 
fied,  a  considerable  number  of  factors  must  be  taken  into 
account.  |i)  The  gradient  energy,  which  is  known  to  be  im¬ 
portant  for  modulations  on  the  scale  considered  here,’-®"’ 
leads  to  dependence  of  the  interdiffusivity  on  the  modulation 
repeat  length.  In  the  present  work,  however,  the  repeat 
length  is  fairly  large,  and  the  reported  interdiffusivities 
are  therefore  probably  close  to  the  bulk  value.  (ii|  Coherency 
strains,  which  are  known  to  affect  the  interdiffusion  in  crys¬ 
tals,'^  are  probably  not  present  in  the  as-prepared  films;  they 
may  become  more  significant,  however,  in  the  course  of  the 
interdiffusion. (iii)  The  composition  dependence  of  the  in- 
lerdiffusion  and  the  attendant  asymmetry  in  the  concentra¬ 
tion  profile  must  be  considered.  (iv|  Effects  of  certain  impuri¬ 
ties  cannot  be  ruled  out  entirely  at  this  stage.  ( v)  The  kinetics 
of  structural  relaxation,  and  its  relation  to  the  pre-existing 
defects  must  be  characterized  in  detail.  The  last  two  factors 
cair  be  elucidated  by  variations  in  the  method  of  preparation. 

At  the  present  stage,  it  can  certainly  be  stated  that  a  new 
diffusion  process  in  amorphous  semiconductors  has  been 
idoitified,  which  is  remarkably  fast,  so  that  it  can  be  mea¬ 
sured  easily  without  risk  of  crystallization  of  the  material. 
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THE  SUPEHMOOULUS  EFFECT  ZN  COKPOSZTIONALLT  MODOLATEO  THIN  FILMS 


R.  C.  Caaanarata* 
Division  of  Applied  Sciences 
Harvard  University 
Cambridge,  Massachusetts  02138 


Introduction 

The  discovery  of  the  supemodulus  effect  in  highly  textured  coirpositionally  modulated  metal¬ 
lic  thin  films  was  sude  by  Hilliard  and  coworkers  (1)  at  Morthwestem  University,  and  most  of  the 
subsequent  ej^rlmental  work  has  been  done  there  (2-9) .  Aie  biaxial  moduli  of  Au-Ni  (1) , 

Cu-Pd  (1,2),  Cu-Ni  (2,3),  and  Ag-Pd  (4,5)  films,  as  measured  by  the  bulge  test  method,  all  showed 
a  single  peak  enhancement  of  between  100-300%  for  a  ccoposition  wavelength  around  20  k,  while 
the  biaxial  modulus  of  Cu-Au  (4,5)  showed  no  such  increase.  Baral  et  al.  (6-8)  reported  a 
sisillar  enhancement  in  the  flexural  modulus,  and  a  bimodal  enhancement  in  Young's  modulus  and  the 
shear  modulus  C^g  (with  the  3-axiE  perpendicular  to  the  plane  of  the  film) .  The  enhancement 
reported  for  the  flexural  modulus  of  (>i-Ni  is  shewn  schematically  in  Fig.  1. 

Some  attempts  to  duplicate  the  above  results  have  been  unsuccessful  (9,10),  though 
Testardi  et  al.  (11)  saw  an  enhancement  in  Young's  modulus  for  a  1^  A  wavelength  film  of  Cu-Ni, 
and  Jankowski  and  Tsakalakos  (12)  Mastired  a  bimodal  enhancement  of  Young's  modulus  in  Cu-NiFe 
films.  A  decrease  in  the  C44  shear  modulus  for  Cu-Hb  (13)  and  Mo-Ni  (14)  thin  films  has  been 
reported  by  Kueny  et  al.  at  the  Argonne  National  Laboratories. 

Theoretical  attespts  to  explain  the  supermodulus  effect  have  been  based  primarily  on  either 
Fermi  surface-Brillouin  zone  interactions  or  coherency  strain  phenomena.  This  article  will 
critically  review  the  current  theoretical  and  experimental  understanding  of  the  supezmodulus 
effect,  and  suggest  avenues  for  further  study. 

Theoretical  Models 

The  major  theoretical  attespts  to  explain  the  supemodulus  effect  can  be  classified  as 
models  in  which  the  elastic  constants  are  affected  by  the  presence  of  either:  (a)  a  Brillouin 
zone  created  by  the  coaposltion  modulation  which  interacts  with  the  Fermi  surface,  causing  a 
dramatic  change  in  the  electronic  structure  of  the  film:  or  (b)  the  relatively  large  coherency 
strains  characteristic  of  textured  modulated  films,  which  introduce  )iigher  order  elastic  effects. 
Any  eoeqplete  theory  of  the  supezmodulus  effect  must  be  able  to  account  for  the  magnitude  of  the 
rather  large  enhancements  (or  dehancementa )  of  the  elastic  constants,  and  also  show  why  the 
effect  only  occurs  near  a  cosposition  wavelength  of  20  A.  In  his  recent  review,  Clapp  (15) 
suBMrily  dismissed  the  theory  based  on  the  coherency  strain  effects,  and  concentrated  on  the 
theories  involving  the  Fermi  surface-Brillouin  zone  interaction,  which  he  claimed  were  the  only 
ones  to  stake  a  serious  atteapt  to  etplain  the  supetsx>dulus  effect.  However,  in  light  of  experi- 
smntal  and  theoretical  work  that  has  appeared  since  his  article  was  presented,  there  are  some 
smjor  difficulties  with  the  Fermi  surface-Brillouin  zone  interaction  theories,  and  the  model 
based  on  coherency  strain  effects  anist  be  considered  to  be  on  at  least  an  equal  footing. 

Fermi  Surface-Brillouin  Zone  Interaction  Theories 

The  idea  of  an  interaction  of  the  Fermi  surface  with  a  Brillouin  zone  sumifesting  itself  in 
the  structure  and  properties  of  solids  is  one  that  dates  back  to  the  Jones  theory  of  alloying 
(16) ,  and  is  found  in  various  forms  throughout  solid  state  physics.  In  almost  every  case  the 
effect  was  originally  described  in  terms  of  second  order  perturbation  theory,  and  the  same  has 

*Current  address:  DepartMnt  of  Materials  Science  and  Engineering,  Massachusetts  Institute  of 
Technology,  Cambridge,  Massachusetts  02139. 
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baan  dona  for  tha  aupazaodulus  thaorlas.  Zt  should  ba  pointed  out  that  the  validity  of  second 
order  perturbation  theory  as  it  has  been  applied  in  these  cases,  as  well  as  the  correctness  of 
accepted  ideas  eoaceming  the  Fermi  surface-Brillouin  zone  interaction  (FSBZI)  in  general,  have 
been  called  into  question  (17) . 

The  enhancement  mechanism  of  the  FSBZZ  is  generally  proposed  in  the  following  way.  A  com¬ 
position  sndulation  with  a  wavevector  ^  creates  artificial  Brillouin  zones  at  1/2  ()c±q) , 
where  li  is  a  reciprocal  lattice  vector.  Assuming  that  one  of  these  artificial  zones  touches 
a  part  of  the  Fersd.  surface  that  is  not  in  contact  with  a  Brillouin  zone  in  the  absence  of  a 
modulation,  the  change  in  electronic  band  structure  due  to  this  critical  contact  could  cause  the 
total  energy  of  a  modulated  film  to  decrease  relative  to  a  hcmogeneous  one.  This  energy  sta¬ 
bilization  will  be  lost  if  the  modulated  foil  is  strained  such  that  a  detachsmnt  of  the  Fermi 
surface  frcm  the  artificial  Brillouin  zone  occurs.  This  detachment  will  manifest  itself  in  a 
stiffer  crystal  cospared  to  one  without  the  cce^osition  sodulation. 

Furdes  (18)  was  the  first  to  suggest  ttiat  if  the  part  of  the  Fermi  surface  that  experienced 
the  critical  contact  was  flat,  a  pseudopotential  calculation  of  the  band  structure  energy  would 
predict  an  energy  stabilization.  Henein  (4)  showed  that  if  effects  due  to  screening  are  ta)ien 
into  account,  an  enhancement  in  the  elastic  constants  is  predicted  due  to  singularities  in  the 
electron  dielectric  function.  A  similar  explanation  has  been  used  to  discuss  Kohn  anomalies 
in  phonon  dispersion  eurvss  of  sooie  elements  (19).  Though  Henein  only  considered  the  case  of  a 


Fcxmi  ^here  mediing  tangential  contact  with  an  artificial  Brillouin  zone,  an  enhancement  viould 
be  elected  for  any  Feziai  surface  topology.  Hu  (20)  has  pointed  out  an  enhancement  could  also 
be  caused  by  singularities  in  the  screened  ion-ion  interaction,  and  that  the  closer  the  modula¬ 
tion  is  to  being  sinusoidal,  the  greater  this  contribution  becoeies. 

The  major  success  that  has  been  claimed  for  these  FSBZl  theories  is  that  X^,  the  wavelength 
of  the  composition  modulation  that  causes  critical  contact,  is  the  wavelength  where  the  enhance¬ 
ment  is  found  experimentally.  However,  given  the  error  involved  in  the  experimental  determina¬ 
tion  of  Fermi  surfaces,  this  succeas  stay  be  overstated.  Taking  the  Cu-Ni  system  as  an  example, 
Hcnein  (4)  and  Hu  (20)  calculated  the  wavelength  necessary  for  critical  contact  (which  is  inde¬ 
pendent  of  the  actual  cause  of  the  energy  stabilization),  ntey  each  used  a  different  Fermi 
surface.  Obtained  from  contradictory  experimental  investigations  (21,22).  As  a  result,  they  each 
))ad  critical  contact  occurring  at  very  different  places  on  the  Fermi  surface.  Despite  these 
totally  different  approaches,  they  both  obtained  a  wavelength  of  critical  contact  close  to  20  X. 
This  suggests  that,  given  the  general  uncertainties  in  Fermi  surface  measurements,  together  with 
the  supermodulus  effect  occurring  at  etbout  the  same  wavelength  independent  of  the  alloy,  any 
agreesaent  between  experiment  and  the  FSBZl  theories  concerning  X^  is  probatbly  fortuitous. 

There  spears  to  be  a  difficulty  with  the  FSBZl  theories  as  formulated  by  Purdes  (18)  and 
Henein  (4) ,  idiere  the  interacting  artificial  Brillouin  zone  planes  are  parallel  to  the  original 
(111)  planes,  in  ea^laining  the  enhancement  found  in  the  shear  modulus  of  Cu-Ni.  Consider 

a  cartesian  coordinate  system  with  the  3-axis  perpendicular  to  the  plane  of  the  film  and  the 
1-  and  2-axes  in  the  plane  of  the  film.  If  0  represents  the  excess  energy  density  created  by 
a  FSBZl,  any  enhancement  6Cgg  in  the  elastic  stiffness  Cgg  is  given  by 


where  Cg  is  the  shear  strain  in  the  plane  perpendicular  to  the  3-axis.  This  means  that  Cg 
must  cause  a  detachment  of  the  Fermi  surface  from  the  artificial  Brillouin  zone  for  there  to  be 
an  enhancenmnt.  This  detachment  could  be  caused  by  the  strain  changing  the  volume  of  the  film, 
thereby  changing  the  Fersii  vector,  and/or  by  the  strain  changing  the  wavelength  of  the  modula¬ 
tion,  thereby  changing  the  reciprocal  lattice  vector  of  the  Brillouin  zone  created  by  the  modu¬ 
lation.  Since  no  voltsne  change  acccnpanies  a  s)iear  strain,  and  since  t))e  shear  strain  Cg  only 
involves  the  1-  and  2-  directions,  the  contact  betMen  the  Fermi  surface  and  the  artificial 
Brillouin  zone  rsnudns  unaffected. 

Based  on  the  above  argwent,  it  seems  that  no  enhaneesMnt  would  be  expected  in  the  shear 
modulus  As  mentioned  in  the  introduction,  Baral  measured  an  enhancement  in  Cgg  on  the 

sane  order  as  the  other  elastic  moduli.  Of  course  this  discrepancy  siay  be  due  to  the  sinplicity 
of  the  model,  and  a  more  detailed  treatment  could  possibly  remedy  tlie  situation.  Nevertheless, 
it  is  disturbing  that  inclusion  of  higher  order  affects  in  the  theory  is  necessary  in  order  to 
explain  an  enhancement  in  that  is  not  needed  for  the  other  moduli. 

It  is  possible  to  argue  within  the  framawoA  of  the  FSBZl  theories  that  no  enhancement  «rould 
be  eigpected  for  the  Cir-Au  system  (20,23).  Zt  is  not  easy,  however,  to  explain  the  ^served  de¬ 
crease  in  the  andulus  Observed  in  Cu-Nb  and  Mo-Ni  foils.  If  cne  extends  Purdes*  and 

Henein* s  theories,  an  enhancement  in  this  elastic  constant  would  be  expected  for  systems  that 
have  a  FSBZl  (24).  Since  Hu* a  theory  is  very  sensitive  to  his  dioice  of  parameters,  it  is  dif¬ 
ficult  to  give  a  precise  description  of  the  prediction  ))is  sodel  would  give;  however,  it  appears 
that  a  general  calculation  would  give  both  an  enhancement  and  a  dehancement  for  different  mve- 
length  regiams  (24) .  (This  enhancesient/dehancement  behavior  is  a  feature  comnon  to  all  calcula¬ 
tions  made  by  Hu  (20).)  Kueny  et  al.  (13)  suggested  that  an  alternative  FSBZl  siechanism  may  be 
involved  in  these  syatems,  though  no  detailed  treatment  was  given. 

Zt  should  be  pointed  out  that  unlike  the  foils  that  displayed  enhancements,  in  which  both 
types  of  layers  were  coc^osed  of  f.c.c.  metals,  the  films  that  showed  a  dehancement  were  made  up 
of  alternating  layers  of  b.c.c.  and  f.c.c.  materials.  As  Baral  has  suggested  (6),  the  elec¬ 
tronic  interactions  in  films  vrlth  ianiscible  components  (such  as  Cu-Nb  and  Mo-Ni)  may  be  signifi¬ 
cantly  different  from  the  filM  that  have  displayed  an  increase  in  the  elastic  moduli.  Thus, 
the  simple  FSBZl  pictuA  may  not  be  so^ilsticated  enough  to  account  for  these  differences. 
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FIG.  2  nie  shear  modulus 
G  as  a  function 
of  cccposition 
modulation  wave* 
length,  as 
calculated  by 
Henein  (4) . 


When  Henein  (4)  performed  a  complete  pseudopotential  calculation  for  the  enhancement  of  a 
shear  modulus  G  (in  a  different  coordinate  system  than  was  used  above)  in  order  to  explain  the 
supermodulus  effect  he  observed  in  Ag-Pd  films,  he  stained  a  result  shown  schematically  in 
Fig.  2.  The  behavior  of  the  curve  comes  from  the  second  derivative  (with  respect  to  the  modula¬ 
tion  wavelength)  of  the  Undhard  screening  dielectric  function,  which  contains  a  logarithmic 
singularity.  As  Henein  pointed  out,  a  negative  shear  modulus  as  obtained  in  Fig.  2  is  a  thermo¬ 
dynamic  ispossibility.  He  argued  that  since  the  wavelength  range  where  the  modulus  was  negative 
was  calculated  to  be  2X10"^  nm,  a  value  of  little  physical  significance,  that  portion  of  the 
figure  could  be  ignored.  However,  a  similar  argument  could  also  be  made  to  ignore  the  enhance¬ 
ment  part  of  the  figure.  The  reason  the  wavelength  range  is  so  small,  is  because  the  Lindhard 
dielectric  function  contains  the  very  wea)c  logarithmic  singular  point.  If  screening  dielectric 
functions  with  stronger  singularities  are  used,  the  range  of  wavelengths  that  display  a  negative 
shear  modulus  will  be  large  enough  that  they  can  no  longer  be  dismissed,  and  this  result  must  Ise 
considered  a  fundasiental  difficulty  of  the  theory.  A  negative  shear  sodulus  is  also  a  possibil¬ 
ity  in  ifu's  model. 

Though  the  FSBZI  theories  are  presently  the  siost  popular  explanations  for  the  supermodulus 
effect,  it  still  remains  to  be  shown  that  these  models  are  able  to  account  for  the  rather  large 
enhancements  in  the  various  moduli  that  have  been  ^served.  The  singularities  on  which  these 
theories  are  based  come  from  the  fact  that  the  calculations  are  done  for  0  K,  thereby  assuming 
a  perfectly  sharp  Fermi  surface.  A  crucial  test  of  the  validity  of  the  FSBZI  theories  is  if 
they  still  predict  large  enhancements  when  a  calculation  is  made  for  a  finite  temperature. 

they  do,  however,  there  is  then  another  problem  to  be  resolved.  Clapp  (15)  made  the  ob¬ 
servation  that  the  types  of  model  tiased  on  an  energy  stabilisation  due  to  Fermi  surface-Brillouin 
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zone  interactions  such  as  have  been  discussed  for  the  supermodulus  effect  have  in  the  past  been 
used  to  explain  the  stability  of  certain  structures.  If  it  is  assumed  that  the  FSBZI  interaction 
indeed  causes  large  enough  energy  changes  to  explain  the  supermodulus  effect,  then  it  would  be 
expected  that  a  modulated  film  with  this  large  energy  stabilization  should  be  thermodynamically 
more  stable  than  a  homogeneous  one.  However,  all  the  foils  that  have  been  displayed  an  enhanced 
modulus  interdiffuse  when  subjected  to  low  tesqperature  anneals,  and  convert  into  homogeneous 
coaposition  films.  It  is  possible,  therefore,  to  argue  that  the  FSBZI  theories  are  doomed  from 
the  Starts  cither  they  will  predict  an  effect  too  aaiall  to  account  for  the  enhancements  observed 
in  the  elastic  stoduli,  or  they  will  predict  a  large  enough  stabilizing  effect  that  is  at  odds 
with  the  observation  that  the  films  are  unstable  to  homogenization. 

Coherency  Strain  Model 

The  other  major  theory  that  has  been  proposed  to  explain  the  supermodulus  effect  is  based  on 
changes  in  the  properties  of  compositionally  modulated  thin  films  due  to  coherency  strain  effects 
(CSE).  The  general  idea  is  that  the  large  coherency  strains  present  in  modulated  films  displace 
the  atoms  from  their  equilibrium  positions  to  such  am  extent  that  they  are  no  longer  within  the 
"Bookean''  part  of  atomic  potential.  As  a  result,  higher  order  elastic  constants  are  necessary  to 
describe  the  elastic  behavior  of  the  film.  These  higher  (third)  order  elastic  constants  are 
manifested  as  enhancements  in  the  standard  second  order  elastic  constants. 

Clapp  (15)  has  argued  that  contributions  to  the  elastic  moduli  due  to  higher  elastic  con¬ 
stants  could  not  even  come  close  to  accounting  for  the  experimentally  observed  enhancements  in 
aK>dulated  films.  However,  Jankows)ii  and  TsakalaXos  (J-T)  have  recently  calculated  such  contri¬ 
butions  in  copper,  for  which  the  third  order  elastic  constants  are  Icnown,  and  showed  that  a  large 
increase  is  possible  (25,26).  They  stated  that,  since  the  symmetry  of  a  film  subjected  to  a 
biaxial  stress  state  (due  to  coherency  strains)  is  changed,  any  calculation  of  elastic  moduli 
must  be  done  using  the  elastic  constants  dictated  by  the  symmetry  of  the  stressed  state.  Unli)ce 
calculations  in  cubic  syametry,  which  gave  no  enhancements  (4,25,26),  J-T  obtained  large  changes 
in  the  (100)  biaxial  modulus  in  copper  determined  in  the  tetragonal  sysmetry  of  its  stressed 
state.  For  example,  if  a  copper  layer  was  subjected  to  a  3%  conpressive  biaxial  strain  due  to 
coherency,  the  biaxial  modulus  would  almost  double  (25,26). 

J-T  also  used  the  pseudopotential  approach  of  Thomas  (39)  to  determine  the  third  order 
contributions  to  the  elastic  constants  of  copper.  The  largest  contribution  to  the  elastic  moduli 
due  to  the  enharmonic  nature  of  the  interatomic  potential  comes  from  the  repulsive  ion-core  over¬ 
lap  energy.  Using  a  suitable  expression  for  this  overlap  energy,  such  as  a  Born-Hayer  potential, 
any  enhancements  in  the  elastic  constants  can  be  calculated  from  the  second  derivative  of  the 
overly  energy  with  respect  to  the  appropriate  strains.  J-T  made  such  calculation  and  obtained 
enhancements  in  the  (100)  biaxial  modulus  similar  to  those  obtained  from  the  steasured  third  order 
elastic  constants  when  oonputed  for  tetragonal  eysanetry.  Calculations  using  the  Born-Mayer 
potentials  for  gold  and  silver  also  showed  large  enhancesienta . 

Because  the  sign  of  the  contribution  from  the  CSE  depends  on  whether  the  stresses  in  the 
layers  due  to  the  coherency  strains  are  compressive  or  tensile  (since  the  contribution  cooies  from 
third  order  terms) ,  such  a  theory  can  potentially  ei^lain  both  the  increases  and  decreases  in 
elastic  snduli  that  have  been  seen  experiaientally.  J-T  also  discussed  that  despite  the  partial 
cancellation  due  to  the  opposite  sign  contributions,  there  will  still  be  a  large  enough  net 
effect  to  explain  the  observed  magnitude  of  the  supermodulus  effect.  However,  it  appears  that 
there  stay  be  some  difficulty  in  explaining  why  no  enhancement  occurs  in  the  cu-Au  system. 

The  major  difficulty  with  the  CSE  theory  is -in  explaining  why  the  enhancements  are  found 
only  for  an  intermediate  range  of  modulation  wavelengths,  or  more  precisely,  why  the  effect  is 
suppressed  for  large  and  small  wavelengths.  J-T  pointed  out  that  for  long  aodulation  wave¬ 
lengths,  it  is  energetically  more  favorable  for  the  film  to  have  adsfit  dislocations  at  the 
interfaces  than  to  retain  the  strain  energy  associated  with  coherency.  If  the  film  does  become 
incoherent  at  long  wavelengths,  it  will  lose  any  effects  due  to  coherency,  such  as  enhancements 
in  the  ixxluli.  Such  an  argument  has  been  used  to  interpret  anomalous  diffusion  behavior  in 
Ctt-Pd  (27).  However,  Oyorgy  et  al.  (28)  have  Inferred  from  their  magnetic  measurements  that 
Cu-Mi  remains  coherent  for  wavelengths  between  16  and  120  A. 


The  arguments  of  J-T  concerning  the  lac)(  of  a  supermodulus  effect  at  small  wavelengths  are 
not  very  convincing.  They  state  that  because  of  ubiquitous  interdiffusion  during  fabrication, 


snail  wavelength  films  do  not  have  as  large  a  ccmposition  modulation  amplitude  as  longer  wave¬ 
length  filns  will.  However,  even  if  the  experimental  measurements  are  corrected  to  take  into 
account  any  differences  in  amplitude,  there  is  still  no  evidence  of  any  enhancement,  as  was  in 
fact  shown  by  "/.akalakos  for  Cu-Mi  (2)  and  J-T  for  Cu-HiFe  (12) .  The  other  argusient  concerning 
small  wavelength  films  offered  by  J-T  (26)  is  not  very  clear:  "...  the  effective  layer  thickness 
for  coherency  strains  across  each  interface  is  about  three  to  five  atomic  dimensions.  Thus  the 
maximum  value  of  the  biaxial  sodulus  occurs  in  a  0.8-1. 2  nm  layer  thiciuiess  (their  italics).” 

This  seems  to  iaply  that  siodulated  films  are  coherent  only  within  a  certain  interfacial  region, 
while  the  rest  of  the  film  has  an  intraplanar  lattice  parameter  that  can  vary  along  the  modula¬ 
tion  direction  (presumably  without  the  introduction  of  dislocations,  which  J-T  have  reserved  for 
longer  wavelength  films).  Allowing  films  to  vary  their  intraplanar  lattice  parameter  in  this  way 
%M3uld  create  astronomically  large  Bhe2Lr  strains,  and  cannot  be  accepted  as  a  realistic  situation. 

TWO  other  points  concerning  str  ill  wavelength  films  are  perhaps  worth  mentioning.  Since 
small  wavelength  (on  the  order  of  two  atomic  spacings)  films  mimic  ordered  structures  while  longer 
wavelength  films  Simulate  phase  separated  structures,  properties  of  modulated  thin  films,  such  as 
the  elastic  moduli,  siay  reflect  these  differences.  Also,  an  elastic  "gradient”  energy  term  has 
been  predicted  to  be  an  is^rtant  contribution  for  small  wavelength  foils  (29)  ,  and  may  affect 
the  effective  elastic  sioduli  measured  in  these  films. 

In  summary,  at  its  current  stage  of  development,  the  CSE  theory  seeam  to  be  able  to  account 
for  the  magnitude  of  both  the  Increases  and  decreases  in  the  elastic  nuaduli  that  have  been  experi¬ 
mentally  determined.  It  is  also  consistent  with  the  observation  that  the  modulated  films  with 
the  supermodulus  effect  are  not  stable  with  respect  to  films  that  have  a  homogeneous  composition. 
Nhat  has  yet  to  be  shown  is  that  the  theory  can  adequately  explain  why  the  supermodulus  effect  is 
only  seen  for  a  limited  range  of  composition  modulation  wavelength.  Despite  this  limitation,  how¬ 
ever,  it  cannot  be  considered  at  this  time  any  less  valid  than  the  FSBZI  theories. 

Since  all  tlie  theories  of  supermodulus  thus  far  presented  are  based  on  quantum  mechanical 
calculations,  it  will  be  necessary  to  consider  the  problem  beyond  the  simple  pseudopotential 
approach  that  has  been  used.  For  example,  it  has  been  shown  that  coherency  strains  have  a  very 
large  effect  on  the  density  of  states  near  the  Fermi  level  (30) .  Though  efforts  are  underway  to 
develop  a  quantum  Mchmically  more  sophisticated  formalism  for  band  structure  determination  in 
modulated  films  (31,32),  no  attempts  have  been  made  to  use  them  to  calculate  enhanced  moduli. 
Obviously,  a  lot  more  work  is  necessary  before  it  can  be  claimed  that  even  a  partial  theoretical 
understanding  of  the  supermodulus  effect  has  been  attained. 

Experimental  Aspects  of  the  Supermodulus  Effect 

Several  methods  have  been  developed  to  measure  various  elastic  smduli  of  modulated  films 
(1,2,6,9,13).  Though  certain  problems  may  exist  with  sene  of  these  tests,  as  will  be  discussed 
below,  it  nevertheless  would  be  desirable  to  employ  all  of  them  on  the  systems  discussed  in  the 
introduction,  since  such  a  catalogue  would  give  a  clearer  empirical  understanding  of  the  super- 
modulus  affect.  Though  it  is  much  sere  easily  said  than  done,  it  would  also  be  worthwhile  to 
develop  a  procedure  to  measure  Poisson's  ratio  V  in  the  plane  of  the  film.  This  modulus  must 
exhibit  a  rather  unusual  dependence  on  the  wavelength  in  Cu-Ni  in  order  to  explain  the  bimodal 
enhancement  of  Young's  modulus  E  and  the  single  peak  enhancements  of  the  bieucial  {e/(1-v)} 
and  the  flexural  {e/(1-v)^}  sK>duli  (2, 3, 6-8).  Baral  (6)  attempted  to  measure  Poisson's  ratio 
in  the  thickness  direction,  and  though  his  results  were  Inconclusive,  they  suggested  that  this 
sndulua  was  negative  for  some  composition  wavelengths. 

Itoxakl  (9)  has  pointed  out  certain  difficulties  associated  with  the  bulge  test,  and  claimed 
that  appto’ent  enhancesmnts  observed  in  the  biaxial  modulus  of  various  nodulated  films  were  just 
experimental  artifacts.  Baral  (6)  has  shown,  however,  that  despite  these  objections  the  super- 
modulus  results  were  still  valid.  Cne  problem  with  the  bulge  test  that  was  not  addressed  was  that 
nonlinear  effects  may  becosm  isportant  for  very  small  asounts  of  deflection,  given  the  thinness 
of  the  films.  This  possibility  has  not  been  incorporated  in  the  analysis  enployed  to  interpret 
bulge  test  results.  A  similar  objection  can  be  raised  with  the  vibrating  reed  method  used  to 
measure  the  flexural  so 'ulus  (24).  Though  this  question  is  important  for  precise  modulus  deter- 
sdnation,  it  probably  does  not  invalidate  the  general  enhanced  modulus  effects  that  have  been 
observed,  since  the  supermodulus  effect  has  demonstrated  a  wavelength  dependence  that  could  not  be 
explained  by  systesiatic  errors  introduced  by  any  nonlinear  aspects  of  the  sothod  of  measurement. 


A  inore  important  issue  concerns  what  modulus  is  actually  measured  in  the  vibrating  reed 
test  (24) .  This  method  involves  measuring  the  resonant  frequency  of  a  film  shaped  as  a  reed, 
clamped  at  one  end.  The  modulus  corresponding  to  this  configuration  depends  on  the  aspect  ratio 
of  the  reed  (width  over  length).  As  this  ratio  goes  to  zero,  young's  modulus  is  measured;  as  the 
ratio  goes  to  infinity,  the  flexural  modulus  is  measured.  For  finite  size  reeds,  St.  Venant's 
principle  suggests  that  the  transition  from  the  one  modulus  to  the  other  occurs  vhen  the  aspect 
ratio  is  close  to  unity,  though  this  would  be  valid  only  for  isotropic  reeds,  tAich  modulated 
filsm  obviously  are  not.  The  foils  measured  by  Baral  were  generally  five  times  longer  than  they 
were  wide  (33),  a  ratio  that  suggests  that  he  actually  measured  Young's  modulus  rather  than  the 
flexural  modulus  as  he  reported.  If  this  is  true,  there  is  a  smjor  discrepancy  laetween  the 
results  for  Cu-Mi  obtained  by  Baral  in  the  vibrating  reed  configuration  and  those  he  reported  as 
Young's  sndulus,  measured  in  a  microtensile  apparatus. 

One  area  of  experimental  investigation  that  requires  more  attention  than  it  has  received  is 
the  relationship  between  the  supermodulus  effect  and  the  structural  nature  of  the  modulated  filnis 
For  example,  it  would  be  interesting  to  see  some  electron  microscopy  done  on  the  interfaces  in 
modulated  films,  such  as  Cu-Ni,  of  various  wavelengths  in  order  to  directly  observe  the  possible 
existence  of  misfit  dislocations.  Such  a  grid  of  dislocations  would  signify  a  loss  of  coherency 
that  is  predicted  to  occur  in  long  wavelength  films,  and  might  explain  changes  in  moduli,  as  was 
discussed  before.  Structural  studies  on  the  systems  that  have  shown  a  modulus  dehancement  are 
also  worth  pursuing.  It  was  reported  (14)  that  the  Mo-Ni  filam  in  tdiich  a  decrease  in  shear 
axjdulus  was  observed  also  showed  evidence  of  a  crystalline  to  amorphous  transition  as  a  fiuction 
of  modulation  wavelength.  This  transition  should  be  characterized  in  more  detail,  since  it 
gives  an  obvious  reason  for  the  abrupt  decrease  in  the  shear  modulus.  Also,  the  35%  decrease 
ih  the  shear  modulus  reported  for  Cu-Nb  films  (13)  (no  precise  value  was  given  for  the  decrease 
in  the  Mo-Ni  foils)  is  characteristic  of  a  crystalline  to  aswrphous  phase  change  in  metals  (34) . 

It  has  been  observed  that  plastic  properties  of  artificially  modulated  structures  have  shown 
a  correlation  with  enhanced  moduli.  For  example,  Baral  et  al.  (6-6)  found  an  enhancement  in  the 
microhardness  of  <hJ-Ni  modulated  films  that  correlated  with  their  vibrating  peed  results.  To 
some  extent  this  is  to  be  expected,  since  elementary  dislocation  theory  predicts  the  energy  of 
edge  and  screw  dislocations  to  be  proportional  to  the  biaxial  and  Young's  moduli,  respectively. 
Isproved  plastic  properties  of  modulated  thin  films  in  fact  were  predicted  by  Koehler  (35)  before 
the  discovery  of  the  supermodulus  effect,  t>ased  on  models  that  considered  only  impediments  to 
dislocation  creation  and  motion,  and  not  on  any  mechanism  to  enhance  the  elastic  B»duli. 

Bevk  et  al.  (36-36)  studied  the  elastic  and  plastic  properties  of  composites  made  up  of  Nb 
filaments  with  diameters  between  SO  and  3000  X  in  a  Cu  matrix.  These  composites  would  be  predic¬ 
ted  to  display  enhanced  resistance  to  plastic  flow,  according  to  Koehler's  theory.  Along  with  an 
increase  in  the  yield  and  tensile  strengths,  an  enhancement  in  the  Young's  modulus  by  a  factor  of 
two  to  three  was  observed.  This  is  quite  intriguing,  since  this  behavior  is  essentially  iden¬ 
tical  with  that  seen  in  Cu-Ni  loodulated  films.  As  mentioned  a)bove,  Koehler's  theory  did  not  in¬ 
volve  any  change  in  the  elastic  moduli,  and  the  theories  of  supermodulus  effect  do  not  seem  to 
apply,  since  the  widths  of  the  filaments  are  oiuch  larger  than  the  modulation  wavelengths  tdiere 
the  enhanced  moduli  have  been  observed  in  modulated  films.  Cohen  and  Bev)(  (37)  suggested  that 
the  enhancement  in  their  composites  nay  originate  from  thezully  induced  elastic  strains  at  the 
Mtrix- filament  interface.  Further  experimental  work  is  warranted  concerning  the  existence  of 
enhanced  moduli  in  composites,  the  effect  being  interesting  in  its  own  right,  but  also  because 
it  might  shed  more  light  on  the  nature  of  the  si^ermodulus  effect  in  compositionally  modulated 
thin  films. 


Ac)cnowiedgement  s 

nie  author  thanks  D.  Baral  for  useful  discussions,  and  F.  Spaepen  and  A.L.  Greer  for  a 
critical  reading  of  the  manuscript.  This  work  was  supported  by  the  Office  of  Naval  Research 
under  contract  N00014-85-K-0023  and  by  the  National  Science  Foundation  under  contract 
DKR83-16979. 


References 


1 

2 

3 


W.M.C.  Yang,  T.  Tsakalakos,  and  J.E.  Hilliard,  J.  Appl.  Phys.  48,  876  (1977). 
T.  Tsakalakos,  Ph.D.  Thesis,  Northwestern  University  (1977). 

T.  Tsakalakos  and  J.E.  Hilliard,  J.  Appl.  Phys.  54,  734  (1963). 


G.  Henein,  t>h.D.  Thesis,  Northwestern  Uhiversity  (1979). 

G.  Henein  end  J.E.  HilUard,  J.  Appl.  Phys.  72B  (19S3). 

D.  Baral,  Ph.D.  Thesis,  Northwestern  University  (1983). 

D.  Baral,  J.B.  Ketterson,  and  J.E.  Hilliard,  in  Modulated  Structure  Materials,  NATO  ASl 
Series,  implied  Sciences,  ed.  T.  Tsa]iala)ios ,  Hartinus  Nijhoff,  Dordrecht,  p.  465  (1984). 

D.  Baral,  J.B.  Ketterson,  and  J.E.  Hilliard,  J.  Appl.  Phys.  1076  (1985). 

H.  Ztozaici,  Ph.D.  Thesis,  Northwestern  university  (1982). 

B.S.  Berry  and  H.  C.  Pritchett,  Thin  Solid  PilsM  33,  19  (1976). 

I. .R.  Testardi,  B.H.  Hillens,  J.T.  Krause,  D.B.  HcWhan,  and  S.  Nakahara,  J.  Appl.  Phys. 

510  (1981). 

A.F.  Jan)iaws)ci  and  T.  Tsakalakoa,  in  layered  Structures,  ^itai^,  and  Interfaces,  Materials 
Research  Society  Sya^sia  Proceedings,  ed.  J.M.  Gibson  and  I>. R.  Dawson,  Materials  Reseaurch 
Society,  Pittsburgh,  Pennsylvania,  Vol.  37,  p.  529  (1985). 

A.  Kueny,  H.  Grinsdltch,  K.  Mlyano,  I.  Banerjee,  C.M.  Falco,  and  Z.X.  Schuller,  Phys.  Rev. 
Lett.  166  (1982). 

M.R.  Khan,  C.S.L.  Otun,  G.P.  Felcher,  M.  Grinsditch,  A.  Kueny,  C.M.  Falco,  and  I.K. 
Schuller,  Phys.  Rev.  B27,  7186  (1982). 

P.C.  Clapp,  in  ref.  7,  p.  455. 

H.  Jones,  Proc.  Boy.  Soc.  A144,  225  (1934). 

V.  Heine  and  D.  Weaire,  in  Solid  State  Physics,  ed.  H.  EhrenreiiA,  F.  Seitz,  and  D. 
Turnbull,  Academic  Press,  NY,  Vol.  24,  p.  250  (1970). 

A.  Purdes,  Ph.D.  Thesis,  Northwestern  University  (1976). 

H.A.  Harrison,  Pseudopotentials  in  the  Theory  of  Metals,  W.A.  Benjamin,  Reading,  HA  (1966). 
T.-B.  Wu,  J.  Appl.  Phys.  53,  5265  (1982). 

M.  Hasegawa,  T.  Suzuki,  and  M.  Hiribayashi,  J.  Phys.  Soc.  Japan  37 ,  85  (1974). 

S.  Nanao,  K.  Kuribayashi,  S.  Taniwaga,  and  M.  Dgyamam,  Phys.  Lett.  38A,  489  (1972). 

W. E.  Pickett,  J.  Phys.  Fs  Met.  Phys.  2195  (1982). 

R. C.  Cammarata,  Ph.D.  Thesis,  Harvard  University  (1985). 

T.  Thakalakos  and  A.F.  Jankowski,  in  ref.  7,  p.  387. 

A.F.  Jankaw8]ci  and  T.  Tsakalakos,  J.  Phys.  F>  Met.  Phys.  15,  1279  (1985). 

A.  Phllofsky,  Ph.D.  Thesis,  Northwestern  University  (1968). 

E. M.  Gyorgy,  J.F.  Dillon,  D.B.  McHhan,  L.H.  tupp,  L.R.  Testardi,  and  P.J.  Flanders,  Phys. 
Rev.  Istt.  45,  57  (1980). 

H.E.  Cook  and  D.  DeFontaine,  Acta  Retail,  607  (1971). 

A.J.  Freeman,  J.  Xu,  and  T.  Jarlborg,  J.  Mag.  6  Mag.  Mat.  31-34,  909  (1983). 

A.  Ctonis  and  N.X.  Flevaris,  Phys.  Rev.  B2S,  7544  (1982). 

A. L.  Yeyatl,  N.V.  Oahan,  and  M.  Weissmann,  Phys.  Rev.  B31,  873  (1985) . 

O.  Baral,  private  coomuni cation. 

D.  Weaire,  M.F.  Ashby,  J.  Icgan,  and  M.J.  Neins,  Acta  Retail.  779  (1971). 

J. S.  Koehler,  Phys.  Rev.  B^,  547  (1970). 

J.  Bevk,  J.p.  Harbison,  and  J.L.  Bell,  J.  Appl.  Phys.  49,  6031  (1978). 

D.E.  Oohen  and  J.  Bevk,  Appl.  Phys.  Lett.  39_,  595  (1981). 

J.  Bevk,  Ann.  Rev.  Mater.  Sci.  319  (1983). 

J.F.  Thomas,  Jr.,  Phys.  Rev.  B^,  2385  (1973). 


39  ,  595  (1981). 


L .-. 

I- •■ 


BASIC  DISTRIBUTION  LIST 


Technical  and  Sutanary  Report? 


Decetober  1982 


Organization 


Codes  Organization 


Copies 


Defense  Docunentation  Center 
Cameron  Station 
Alexandria,  VA  22314 

Office  of  Naval  Research 
Department  of  the  Navy 
800  N.  Quincy  Street 
Arlington,  VA  22217 
Attn;  Code  431 

Naval  Research  Laboratory 


12 


Naval  Air  Propulsion  Test  Center 
Trenton,  NJ  08628 
ATTN:  Library 

Naval  Construction  Battallion 
Civil  Engineering  Laboratory 
Port  Hueneme,  CA  93043 
ATTN:  Materials  Division 


Naval  Electronics  Laboratory 


Washington,  DC  20375 

San  Diego,  CA  92152 

I 

ATTN;  Codes  6000 

1 

ATTN:  Electron  Materials 

6300 

1 

Sciences  Division 

-V 

U*  • 

2627 

1 

>  , 

Naval  Air  Development  Center 

Naval  Missile  Center 

Code  606 

Warminster,  PA  18974 
ATTN;  Dr.  J.  DeLuccia 

Commanding  Officer 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 
ATTN:  Library 

Naval  Oceans  Systems  Center 
San  Diego,  CA  92132 
ATTN;  Library 

Naval  Postgraduate  School 
Monterey,  CA  93940 
ATTN;  Mechanical  Engineering 
Department 

Naval  Air  Systems  Command 
Washington,  DC  20360 
ATTN;  Code  31A 

Code  5304B 


Naval  Sea  System  Command 
Washington,  DC  20362 
ATTN:  Code  05R 


Materials  Consultant 

Code  3312-1 

Point  Mugu,  CA  92041 

Commander 

David  W.  Taylor  Naval  Ship 
Research  and  Development  Center 
Bethesda,  MD  20084 


Naval  Underwater  System  Center 
Newport,  RI  02840 
ATTN:  Library 

Naval  Weapons  Center 
China  Lake,  CA  93555 
ATTN:  Library 


NASA 

Lewis  Research  Center 

21000  Brookpark  Road 

Cleveland,  OH  44135 

ATTN;  Library  1 

National  Bureau  of  Standards 

Washington,  DC  20234 

ATTN:  Metals  Science  and  Standards 

Division  1 

Ceramics  Glass  and  Solid  State 
Science  Division  1 

Fracture  and  Deformation  Div.  1 


*  .  ^ 

■ ,  ^ 


L-r- 


•s*  ■ 


RE/431/82/174 


I 


Naval  Facilities  Engineering 
Command 

Alexandria,  VA  22331 

ATTN:  Code  03  1 


Scientific  Advisor 

Commandant  of  the  Marine  Corps 

Washington,  DC  20380 

ATTN:  Code  AX  1 

Army  Research  Office 
P.  0.  Box  12211 
Triangle  Park,  NC  27709 
ATTN:  Metallurgy  &  Ceramics 

■ '  Program  1 


Army  Materials  and  Mechanics 
Research  Center 
Watertown,  MA  02172 
ATTN:  Research  Programs 
Office 

Air  Force  Office  of  Scientific 
Research/NE 
Building  410 
Bolling  Air  Force  Base 
Washington,  DC  20332 
ATTN:  Electronics  &  Materials 
Science  Directorate  1 


Defense  Metals  and  Ceramics 
Information  Center 
Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  OH  43201 

Metals  and  Ceramics  Division 
Oak  Ridge  National  Laboratory 
P.O.  Box  X 

Oak  Ridge,  TN  37380 

Los  Alamos  Scientific  Laboratory 
P.O.  Box  1663 
Los  Alamos,  NM  87544 
ATTN:  Report  Librarian 


Argor.ne  National  Laboratory 
Metallurgy  Division 
P.O.  Box  229 
Lemont,  IL  60439 


Brookhaven  National  Laborato'"' 
Technical  Information  Division 
Upton,  Long  Island 
New  1£ork  11973 
ATTN:  Research  Library 


Library 

Building  50,  Room  134 
Lawrence  Radiation  Laboratory 
Berkeley,  CA 


NASA  Headquarters 

Washington,  DC  20546 

ATTN:  Code  RRM  1 


Send  ONE  COPY  to  each  unless  otherwise  indicated. 


